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Abstract 

Hadronic event shape distributions from e + e~ annihilation measured by the OPAL ex- 
periment at centre-of-mass energies between 91 GeV and 209 GeV are used to determine 
the strong coupling ag. The results are based on QCD predictions complete to the next- 
to-next-to- leading order (NNLO), and on NNLO calculations matched to the resummed 
next-to- leading- log- approximation terms (NNLO+NLLA). The combined NNLO result 
from all variables and centre-of-mass energies is 

a s (m z0 ) = 0.1201 ± 0.0008(stat.) ± 0.0013(exp.) ± 0.0010(had.) ± 0.0024(theo.). 

while the combined NNLO+NLLA result is 

a s (m z0 ) = 0.1189 ± 0.0008(stat.) ± 0.0016(exp.) ± 0.0010(had.) ± 0.0036(theo.). 

The completeness of the NNLO and NNLO+NLLA results with respect to missing higher 
order contributions, studied by varying the renormalization scale, is improved compared 
to previous results based on NLO or NLO+NLLA predictions only. The observed energy 
dependence of ag agrees with the QCD prediction of asymptotic freedom and excludes 
the absence of running. 
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1 Introduction 



Events originating from e + e~ annihilation into quark- ant iquark pairs allow precision 
tests [TH3] of the theory of the strong interaction, Quantum Chromodynamics (QCD) jl]- 
[7J. Comparison of observables like jet production rates or event shape variables with 
theoretical predictions allows the crucial free parameter of QCD - the strong coupling ag 
- to be determined. The determination of as from many different observables provides an 
important consistency test of QCD. Recently, significant progress in the theoretical calcu- 
lation of event shape observables has been made. Next-to-next-to-leading order (NNLO) 
calculations are now available jHllS] as well as NNLO calculations matched with resummed 
terms in the next-to- leading- logarithmic-approximation (NLLA) |lUj . 

Measurements of as at centre-of-mass-system (cm.) energies between y/s = 91 GeV 
and 206 GeV, using NNLO predictions and ALEPH data, were presented in [11]. This 
was followed by measurements between y/s = 14 GeV and y/s = 44 GeV [12] and between 
y/s = 91 GeV and y/s = 206 GeV [T3] based on comparing JADE or ALEPH data with 
NNLO plus matched NLLA calculations. The strong coupling has also been determined 
at NNLO from the three-jet rate at LEP [14] . In the present study, the revised NNLO 
calculations described in [8l|9] are used to determine as at NNLO and NNLO+NLLA for 
13 different energy values between 91 and 209 GeV. The data sample is that of the OPAL 
Collaboration at LEP. Hadronization corrections are treated in the manner described 
in [2JIHIHZIHS] - The study is based on the same measurements of event shape variables, 
the same fit ranges (except for yf 3 as discussed below), and the same Monte Carlo models 
for detector and hadronization corrections as those in [15]. 

The structure of the paper is as follows. In Sect. [2] we give an overview of the OPAL 
detector, and in Sect. |3] we summarize the data and Monte Carlo samples used. The 
theoretical background to the work is outlined in Sect. HI The experimental analysis 
techniques are explained in Sect. |5j and the measurements are compared with theory in 
Sect. El 



2 The OPAL Detector 

The OPAL experiment operated from 1989 to 2000 at the LEP e+e~ collider at CERN. 
The OPAL detector is described in detail in [ToT - fTS] . This analysis mostly makes use of 
the measurements of energy deposited in the electromagnetic calorimeter and of charged 
particle momenta in the tracking chambers. 

All tracking systems were located inside a solenoidal magnet, which provided a uniform 
axial magnetic field of 0.435 T along the beam axis. The main tracking detector was the 
central jet chamber. This device was approximately 4 m long and had an outer radius of 
about 1.85 m. The magnet was surrounded by a lead glass electromagnetic calorimeter 
and a sampling hadron calorimeter. The electromagnetic calorimeter consisted of 11704 
lead glass blocks, covering 98% of the solid angle. Outside the hadron calorimeter, the 
detector was surrounded by a system of muon chambers. 



5 



3 Data and Monte Carlo Samples 



The analysis is based on the same data, Monte Carlo samples, and event selection, as 
those employed in a previous study [T5]. The data were collected at cm. energies be- 
tween 91.0 GeV and 208.9 GeV. The data at 91 GeV, on the Z° peak, were collected 
during calibration runs between 1996 and 2000 and have the same conditions of detector 
configuration, performance and software reconstruction as the higher energy data. The 
data are grouped into samples with similar cm. energies, as indicated in Tab. [IJ Be- 
sides the energy range, mean energy and integrated luminosity, Tab. [T] lists the year(s) of 
collection and the number of selected hadronic annihilation events for each sample. 

Samples of Monte Carlo-simulated events were used to correct the data for experi- 
mental acceptance, efficiency and backgrounds. The e + e~ —> qq process was simulated 
at V« = 91 - 2 GeV usin S JETSET 7.4 [T3], and at higher energies using KX2f 4.01 or 
/C/C2f 4.13 [2D1EU with hadronization performed using PYTHIA 6.150 or PY- 
THIA 6.158 pE]. Corresponding samples of HERWIG 6.2 [221123] or /COf events with 
HERWIG 6.2 hadronization were used for systematic checks. Four-fermion background 
processes were simulated using grc4f 2.1 [24], KORALW 1.42 [23] with grc4f [21] matrix 
elements, with hadronization performed using PYTHIA. The above samples, generated 
at each energy point studied, were processed through a full simulation of the OPAL 
detector [26] and reconstructed in the same manner as the data. In addition, when 
correcting for the effects of hadronization, large samples of generator-level Monte Carlo 
events were employed, using the parton shower models PYTHIA 6.158, HERWIG 6.2 and 
ARIADNE 4.11 [27]. 

Each of these models used for the description of hadronization and detector response 
contains a number of tunable parameters. These parameters were adjusted to describe 
previously published OPAL data at ^Js ~ 91 GeV as discussed in [2B] for PYTHIA/ 
JETSET and in [29] for HERWIG and ARIADNE. 



4 Theoretical Background 
4.1 Event Shape Distributions 

The properties of hadronic events can be described by event shape observables. The 
event shape observables used for this analysis are thrust (1 — T) [301131], heavy jet mass 
(Mh) [32], wide and total jet broadening (5w and B^) [33], C-Parameter (C) [34H36] 
and the transition value between 2 and 3 jet configurations defined using the Durham jet 
algorithm (y^ 3 ) [3?]. 

In the following, the symbol y is used to refer to any of the variables 1 — T, Mh, -Bt, 
£>w ; C or y® 3 . Event shape variables characterize the main features of the distribution of 
momentum in an event. Larger values of y correspond to the multi-jet region dominated 
by the radiation of hard gluons. Smaller values of y correspond to the two-jet region with 
only soft and collinear radiation. 
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4.2 QCD Calculations 

QCD predictions for the distribution of event shape observables in e + e~ annihilations are 
now available to 0(a^) (NNLO) [SUH]- In the case where the renormalization scale [Ir 
equals the physical scale Q identified with ^/s, the predictions are: 



1 da 
cr dy 



0(a|) 



dA dB 2 dC , 
d^ S + d^ S + crT s 



with as = as(nR)/(27r). The coefficient distributions for the leading order (LO) dA/dy, 
next-to-leading order (NLO) dB/dy and NNLO dC/dy terms were provided to us by 
the authors of [8]. The normalization to the total hadronic cross section and the terms 
generated by variation of the renormalization scale parameter = Hr/Q are implemented 
according to the prescription of [8]. 

In the two-jet (low y) region, the effect of soft and collinear emissions introduces large 
logarithmic effects depending on L = log(l/y). For the (1 — T), M H , B T , B w , C and 
y® 3 distributions, the leading and next-to-leading logarithmic terms can be resummed up 
to infinite order in perturbation theory [3B]. This is referred to as the next-to-leading- 
logarithmic approximation (NLLA). The most complete calculations of event shape ob- 
servables are obtained from combining (matching) the O(o,g) and NLLA calculations, 
taking care not to double count terms that are in common between them. 

The matching is not unique. In the so-called InR matching scheme, the NNLO+NLLA 
expression for the logarithm of the cumulative distribution R(y) = f dy'(l/a)da/dy' 
is [10] 

In (R(y,6i s )) = Lgx (a s L) + g 2 (6i s L) (2) 
+ a s (A(y)-GnL-G 12 L 2 ) 

+ al^B(y)-^A 2 (y)-G 22 L 2 -G 23 L 3 
+ &l (c (y) — A(y)B (y) + ^ 3 (y) - G 33 L 3 - G M V 

The functions g\ and g 2 represent the resummed leading and next-to-leading logarithmic 
terms while the Gij are matching coefficients. The coefficient functions A, B and C 
are related to the differential coefficients in (TT]) by integration, e.g. A(y) = dy' d - . 
Since as in the NNLO and NLLA terms are assumed to be the same, there is only one 
renormalization scale. The NLLA terms introduce a further arbitrariness in the choice of 
a logarithmic rescaling variable x L , see [T5j . 

The theoretical calculations provide distributions at the level of quarks and gluons, 
the so-called parton-level. Monte Carlo-based distributions calculated using the final- 
state partons after termination of the parton showering in the models are also said to be 
at the parton-level. In contrast, the data are corrected to the hadron-level, i.e. they corre- 
spond to the distributions of the stable particles in the event as explained in Sect. 15.21 To 
compare the QCD predictions with measured event shape distributions these predictions 
are corrected from the parton to the hadron level. The corrections are based on large 
samples of typically 10 7 events generated with the parton-shower Monte Carlo programs 
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PYTHIA (used by default), HERWIG and ARIADNE (used to evaluate systematic uncer- 
tainties). The corrections are defined by the ratio of the Monte Carlo distributions at the 
hadron and parton levels and are applied as multiplicative corrections to the theoretical 
predictions. 

We compare the parton-level calculations of the Monte Carlo generators with the 
QCD calculations in NNLO+NLLA with ag(mzo) = 0.118, x M = 1.0 and x L = 1.0 at 
y/s = 91 GeV, see Fig. [TJand [12]. The dashed line shows the deviation from one of the 
ratio between the Pythia prediction at the parton level and the NNLO+NLLA calculation. 
The dotted (dash-dotted) line shows the corresponding result from Herwig (Ariadne). 
The solid lines, symmetric about zero, show the maximum difference between any pair of 
the Monte Carlo models, which are seen to be of similar size to the differences between 
the generators and the NNLO+NNLA calculation. Therefore we consider that the MC 
simulations adequately account for the hadronization correction to the calculations. 

The model dependence of the hadronization correction is included as a systematic 
uncertainty as described below. 



5 Experimental Procedure 

The data used in the present paper are identical to those presented in [15],[39j|l0]. For 
completeness and to facilitate the discussion of systematic uncertainties, we give a brief 
summary of the analysis procedure below. 

5.1 Event Selection 

The event selection procedure is described in [15]. This procedure selects well- measured 
hadronic event candidates, removes events with a large amount of initial-state radiation 
(ISR) at 130 GeV and above, and removes four-fermion background events above the 
W + W~ production threshold of 160 GeV. 

5.2 Corrections to the Data 

For each accepted event, the value of each event shape observable is computed. A standard 
algorithm [JT] is applied to mitigate the double-counting of energy between the tracking 
chambers and calorimeter. To correct for background contributions, the expected number 
of residual four-fermion events hi is subtracted from the number of data events in each 
bin i of each distribution. Simple bin-by-bin corrections are applied to account for detector 
acceptance, resolution, and the effects of residual ISR. We examine the Monte Carlo 
predictions for the event shape variables at two levels: the detector level, which includes 
simulation of the detector and the same analysis procedures as applied to the data, and 
the hadron level, which uses stable particles^ only, assumes perfect reconstruction of the 
particle momenta, and requires the hadronic cm. energy yfs' to satisfy \fs — y/P < 1 GeV. 
The ratio of the hadron- to the detector-level prediction for each bin, ct^, is used as 
the correction factor for the data, yielding the corrected bin content Ni = oii(Ni — hi). 

1 For this purpose, all particles having proper lifetimes greater than 3 x 10 -10 s are regarded as stable. 
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This corrected hadron-level distribution is then normalized by the total number of events 
N = E fc N k and bin width W t : Pi = ^/(JWj). 

5.3 Systematic uncertainties 

To evaluate systematic uncertainties, the analysis is repeated after modifying the selection 
or correction procedures. The event and track selection cuts are varied within the ranges 
of values suggested in [T5]. The Monte Carlo model employed to calculate the detector 
correction is altered, and the cross section used in the subtraction of four-fermion events 
is varied by ±5%. In each case, the difference in each bin with respect to the standard 
analysis is taken as a contribution to the systematic uncertainty 

6 Determinations of as 
6.1 Fit procedure 

The strong coupling a§ is determined using a minimum-x 2 fit, comparing theory with each 
of the measured event shape distributions at the hadron level. A x 2 value is calculated 
at each cm. energy: 



where i, j include the bins within the fit range of the event shape distribution, Pj is the 
measured value in the ith bin, £j(as) is the QCD prediction for the ith bin corrected for 
hadronization effects, and V~ l is the inverse of the statistical covariance matrix V stat of 
the values P;. The QCD prediction is obtained by integrating the result in ([ID over the 
bin width, and then applying the hadronization correction. The \ 2 value is minimized 
with respect to as with the renormalization scale factor x M and the rescaling variable x L 
set to 1. The evolution of the strong coupling as as a function of the renormalization scale 
is implemented to three-loop order |42j . Since the cm. energy range considered here does 
not cross flavour thresholds, no significant uncertainties are introduced by the evolution 
of as- Separate fits are performed to each of the six observables at each cm. energy value. 

To account for correlations between bins in the computation of x 2 ? the statistical 
covariance matrix V^ tat is calculated following the approach described in [T5] : 



where So- is the Kronecker delta function, and the other terms are defined in Sect. 15.21 

To test our procedure, we verified that we are able to reproduce the fit results from [151 
H3] using NLO+NLLA predictions within typically 0.5%. 




(3) 




(4) 



k 
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The fit ranges are chosen as in [T5ll4*4"] and correspond to regions where the corrections 
of the data and of the parton level theory to the hadron level are both reasonably small 
and where the NLO+NLLA fit results are stable under small variations of the fit range. 
The relevant calculations in the two-jet region are the NLLA terms, which are identical 
in the present study and [TS]. The NNLO prediction for y^ 3 is given in bins of — lnyf 3 . 
The transformation to bins of y^ 3 produces precision problems for small values of yf 3 . 
Therefore we use a more restrictive fit range for y^ 3 than in [15], which approximately 
matches the range in [12]: the lower bound is 0.012 rather than 0.0023. 

The fit ranges are shown in Tab. [2j In the resummation of log-enhanced terms, the 
leading log term of dA/dy is \ia(y)/y. At y values equal to the lower limit y , the value 
of a$ \n(y )/y is still of the order of one for as(m z o) = 0.118 and 91 < yfs < 207 GeV. 
Thus the fit ranges are also suitable for pure NNLO analysis. 

Previously measured event shape distributions are published [15] only in the energy 
ranges as separated by horizontal lines in Tab. [TJ We perform fits to extract as at all of 
the energy points^ shown in Tab. [TJ column 2. We average over the larger energy ranges 
only for purposes of presentation. 

The detector correction factors are typically between 0.9 and 1.5 within the fit ranges. 
For energies above 189 GeV, which have low statistics, the corrections exceed 2.0 in the 
multi-jet regions. The maximum correction of 4.8 occurs for the 1 — T distribution at 
207 GeV. The variable with the least variation in the size of its hadronization correction, 
and with hadronization corrections closest to one, is y^ 3 . 

The statistical uncertainty of as is given by the variation required to increase \ 2 by 
one unit from its minimum. The systematic uncertainties account for experimental effects, 
the hadronization correction procedure and uncertainties of the theory. The three sources 
of systematic uncertainty are added in quadrature to the statistical uncertainty to define 
the total uncertainty. Below, we describe the evaluation of systematic uncertainties. For 
each variant of the analysis, the corresponding distribution is fitted to determine as, and 
the difference with respect to the value of as from the default analysis is taken as a 
systematic uncertainty contribution. 

Experimental Uncertainties: These are assessed as described in Sect. 15.31 

Hadronization: For the default analysis, PYTHIA is used to evaluate hadronization 
corrections (Sect. |4T2~|) . As systematic variations, HERWIG and ARIADNE are used 
instead. The larger of the deviations is taken as the systematic uncertainty. It was 
observed in [451 H6] that systematic uncertainties determined from the differences 
between the PYTHIA, HERWIG, and ARIADNE models are generally much larger 
than those that arise from varying the parameters of a given model. 

Theoretical Uncertainties: The theoretical calculation of event shape observables is a 
finite power series in as- The uncertainties originating from missing higher order 
terms are assessed by changing the renormalization scale factor to x M = 0.5 and 
= 2.0. The rescaling variable is set to x L =2/3 and x L = 3/2 (x L =4/9 and 
x L = 9/4 in case of y^ 3 ) [15]. The largest deviation with respect to the standard 

2 This allows the use of the event shape covariance matrices which are energy dependent and have been 
calculated separately for each energy point. 
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analysis is taken as the systematic uncertainty. A variation of the matching scheme 
is not studied because i?-matching is not available for NNLO+NLLA. The match- 
ing scheme was not an important source of uncertainty in studies of as based on 
NLO+NLLA calculations. 

6.2 Results from NNLO Fits 

The results of the NNLO fits are summarised in Tab. [3j In order to clarify the presenta- 
tion, the as values from 130 and 136 GeV are combined^], and likewise the values from 
161 to 183 GeV, and from 189 GeV and above. These three ranges of cm. energies, 
which correspond to luminosity weighted mean energies of 133.1, 177.4 and 197.0 GeV, 
respectively, cover sufficiently small ranges of y/s that only small variations of as occur 
within each range. 

Figs. El [3] and S] show the 1 — T, Mr, Bt, -Bw, C and y® 3 event shape distributions 
together with the NNLO fit results for the 91.3 and 206.6 GeV energy points. All event 
shape distributions are described reasonably well by the fitted predictions. The x 2 /d.o.f. 
values, which are large for the data at the Z° peak because they are based on statis- 
tical uncertainties only, range from 0.03 for y® 3 at \fs = 191.6 GeV to 228 for Bt at 
\fs = 91.3 GeV. The uncertainties at this point are dominated by the experimental sys- 
tematic uncertainties (discussed below). To check how the inclusion of the experimental 
errors in the fits reduces the large \ 2 values, an alternative \ 2 value is defined as fol- 
lows. An estimate of the experimental systematic covariance is added to the statistical 
covariance (JJj), 

^total = yg* + min(a 2 xp i; a 2 xpj); (5) 

where cr exPi j is the experimental systematic uncertainty at bin %. Tab. H] shows the fit 
results at 91 GeV. The fit results are in general not compatible within the combined 
statistical and experimental errors with the results when only the statistical covariance 
is used. These y^/d.o.i. values are smaller than those from employing the statistical 
covariances, however they still range up to 62 for Bt- We return to this issue when we 
discuss the NNLO+NLLA fits below. 

The results for as from the different event shape variables are remarkably consistent 
with each other, for each cm. energy, as seen from Tabs. [TQl [TT] (but not in the case of 
the modified fit procedure whose results are shown in Tab. H]). We find root-mean-square 
(r.m.s.) values for a s (y/s) between 0.0012 at 183 GeV and 0.0044 at 202 GeV. The 
variations in as between the different variables are comparable to the total systematic 
uncertainties. The values of as are significantly larger at 91 GeV than at higher energies, 
providing evidence for the running of as- Theoretical uncertainties dominate at the 
Z° peak (except for y® 3 ), where the statistical uncertainties are small. Similarly, statistical 
uncertainties are small at 189 GeV where there is a relatively large data sample. Statistical 
uncertainties dominate at 130, 161 and 172 GeV, where the data samples are smaller. 
Experimental systematic uncertainties are small at 91 GeV and larger at higher y/s where 

3 This combination is performed assuming minimum overlap correlation of the experimental uncertain- 
ties at the different cm. energy points. The procedure includes a correction for the running of as(\/s) 
to the luminosity weighted mean energy. 
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uncertainties from subtraction of ISR and four fermion events contribute more strongly. 
For completeness, the fit results for all energy points are given separately in App. |Aj 

6.3 Results from NNLO+NLLA Fits 

The results of the NNLO+NLLA fits are shown in Figs. [2, E] and |4] and listed in Tab. [51 
The calculations fit the data better than the pure NNLO calculations at small values of 
the variables. The values of x 2 /d.o.i. are smaller on average than for the NNLO fits, 
especially at 91 GeV, indicating better consistency with the data. The r.m.s. values of 
<*s(y/s) vary between 0.0017 at yfs = 205 GeV and 0.0051 at -y/s = 202 GeV, i.e. the 
scatter of individual results is essentially the same as for the NNLO analysis. The pattern 
of statistical, experimental, and hadronization uncertainties is the same as for the NNLO 
fits discussed above. At most energy points the theory error results from the x L variation 
for Mh, -Bw, y® 3 , and from the x M variation for 1 — T, B^, C. Compared with the NNLO 
analysis the values of as are lower by 0.4% on average, and the theoretical uncertainties 
are higher by 33%. Larger scale uncertainties are expected to arise at 91 GeV when the 
NLLA terms are added to the NNLO terms [10l[T3j|17], because the NNLO calculation 
compensates for the variation of the renormalization scale in two loops, while the NLLA 
term compensates for the variation in only one loop. The difference we observe between 
as in the NNLO and NNLO+NLLA studies is smaller than the corresponding difference 
observed between the NLO and NLO+NLLA studies (Sect. 1675]) . as predicted in [10]. This 
difference is also smaller than that observed at lower energies [12], as expected from the 
energy evolution of as- 

The fit results at 91 GeV employing the total covariance © are shown in Tab. |6j 
Unlike the case of pure NNLO, the fit results are compatible with the results based on 
only the statistical covariance within the combined statistical and experimental systematic 
uncertainties. The x 2 /d.o.i. values are of the order of 1. For a reasonable description of 
the data by the theory predictions, inclusion of the resummed logarithmic terms (NLLA) 
is seen to be important. As the covariance (EJ) is only an approximation, we do not use 
the fit values or errors further. 

6.4 Combination of Results 

The results obtained at each energy point for the six event shape observables are com- 
bined using uncertainty- weighted averaging as in [2"}[T5' } |4"5 } |4"9"] . The statistical correlations 
between the six event shape observables are estimated at each energy point from fits to 
hadron-level distributions derived from 50 statistically-independent Monte Carlo samples. 
The experimental uncertainties are determined assuming that the smaller of a pair of cor- 
related experimental uncertainties gives the size of the fully correlated uncertainty. The 
minimum overlap assumption results in a conservative estimate of the total uncertainty. 
The hadronization and theoretical systematic uncertainties are evaluated by repeating 
the combination with changed input values, i.e. using a different hadronization model 
or the different value of x^ or x L that yields the maximum deviation. The results are 
given in Tab. [7] and shown for the NNLO+NLLA analysis in Fig. [5j Also shown is a 
comparison with the results of a NNLO+NLLA analysis at lower energy from the JADE 
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Collaboration [T2] . 

To study the compatibility of our data with the QCD prediction for the evolution 
of the strong coupling with cm. energy we repeat the combinations with or without 
evolution of the combined results to the common scale, setting the theory uncertainties to 
zero since these uncertainties are highly correlated between energy points. We assume the 
hadronization uncertainties to be partially correlated. The \ 2 probability of the average 
for a running; (constant) coupling in the NNLO+NLLA study is 0.59 (6xl(r 13 ). The 
corresponding result for the NNLO study is 0.56 (2xl0 -6 ). We interpret this as clear 
evidence, from OPAL data alone, for the running of as in the manner predicted by QCD. 

We evolve the results in Tab. [7| with \fs > m%o to a common scale myp using the 
QCD formula for the running of as- The results are then combined using the uncertainty- 
weighted averaging procedure mentioned above. The results are shown in Tab. |H] for 
the NNLO and NNLO+NLLA analyses. The consistency between these values and the 
results from 91 GeV (Tab. [7]) demonstrates the compatibility of the data with the running 
predicted by QCD. The high energy data have smaller theoretical and hadronization 
uncertainties, and therefore complement the statistically superior 91 GeV data. 

The values of as(r«z ) obtained from combining all energy points are also given in 
Tab. E The resulting values of a s (m z o) = 0.1201 ± 0.0031 from the NNLO study, and 
a s (m z o) = 0.1189 ± 0.0041 from the NNLO+NLLA study, are consistent with the world 
average [50] (0.1184 ± 0.0007), recent NNLO analyses of JADE [TJ] (0.1172 ± 0.0051) 
and ALEPH [TJ] (0.1240 ± 0.0033) data, our previous study [15] based on NLO+NLLA 
calculations (0.1191 ± 0.0047) and the study p3] of ALEPH data using NNLO+NLLA 
(0.1224 ± 0.0039). The total uncertainties of 2.6% and 3.4% for ag(m z o) place these 
measurements amongst the most precise determinations of as available. 

After running the fit results for «s(v^) f° r each observable to the common reference 
scale m z o, we combine the results for a given observable to a single value. We use the 
same method as above and obtain the results for as(?n z o) shown in Tab. [91 These results 
demonstrate that the measurements from the different observables are far from compatible 
with each other when only statistical uncertainties are considered, but are consistent with 
a common mean when the systematic uncertainties are included, neglecting their correla- 
tion. The r.m.s. values of the results for as(?«z ) are 0.0013 for the NNLO analysis and 
0.0026 for the NNLO+NLLA analysis; both values lie within the range of the uncertainty 
for the corresponding combined result shown in Tab. [HJ Fig. [6] displays the combined as 
result for each observable. Results from NLO and NLO+NLLA studies, discussed below, 
are also shown. Combining the NNLO or NNLO+NLLA results of Fig. [6] to obtain an 
overall value for as, or evolving each event shape measurement from each energy to the 
reference scale and then combining, yields the same result as given by the corresponding 
measurement in Tab. [8] to within Aas( m z°) = 0.0003. 

6.5 Comparison with NLO and NLO+NLLA fits 

To compare our results with previous as measurements, the fits to the event shape distri- 
butions are repeated with NLO predictions and with NLO predictions combined with 
resummed NLLA with the modified In i?-matching scheme (NLO+NLLA), both with 
i„ = 1. The NLO+NLLA calculations with the modified In R- matching scheme were the 
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standard of the analyses at the time of termination of the LEP experiments [T5 | l4*9 | [5T | l52] . 
The fit ranges and procedures for the evaluation of systematic uncertainties are the same 
as those used above for the NNLO and NNLO+NLLA studies and thus differ somewhat 
from the previously published results [TBI 11)5] . 

The combination of the fits using NLO predictions yields as(^z°) — 0.1261 
±0.0011 (stat.) ± 0.0024(exp.) ± 0.0007(had.) ± 0.0066(theo.) while the combination of 
NLO+NLLA results yields a s (m z o) = 0.1173 ±0.0009(stat.) ±0.0020(exp.) ±0.0008(had.) 
±0.0055(theo.). These results are shown by the corresponding shaded bands in Fig. |6j 
The result obtained with the NLO+NLLA prediction is consistent with the NNLO and 
NNLO+NLLA analyses, but the theory uncertainties are larger by a factor of 2.3 and 1.5 
respectively. The analysis using NLO predictions gives theoretical uncertainties larger by 
a factor of 2.8 and 1.8, and the value for ag(mzo) is larger compared to the NNLO or 
NNLO+NLLA results. It has been observed previously that values for ctg from NLO 
analyses with large in comparison with most other analyses [51]. The 

NLO+NLLA analysis yields a smaller value of as{ m z°) compared to the NLO result, 
and the NNLO+NLLA analysis a smaller one compared to NNLO. The difference be- 
tween the NNLO+NLLA and NNLO results is smaller than the difference between those 
of the NLO+NLLA and NLO fits because a larger part of the NLLA terms is included in 
the NNLO calculations. 

6.6 Renormalization Scale Dependence 

To assess the dependence of as on the choice of the renormalization scale, the fits to 
distributions of the six event shape variables are repeated using NNLO, NNLO+NLLA, 
NLO and NLO+NLLA predictions with 0.1 < x^ < 10. As an example, as and the 
X 2 /d.o.f. for the C variable at \/s = 91 GeV are shown as a function of x^ in Fig. The 
smallest x 2 /d.o.f. values for the fixed order calculations arise at the smallest scales, while 
with the NLLA calculations smaller x 2 /d.o.f. values occur nearer the physical scale y/s, 
i.e. Xfj,—1. Near this scale, smaller x 2 /d.o.f. values are observed with the NNLO curves 
than with the respective NLO curves. 

The NLO calculation yields a larger value of as than the other calculations for x M > 
0.2. The a s (m z o) values using the NLO+NLLA, NNLO and NNLO+NLLA calculations 
cross near the naturafl choice of the renormalization scale x M = 1. The NLLA terms at 
Xfj, — 1 averaged over the fit range are almost identical to the (9(ag)-terms in the NNLO 
calculation. A similar behaviour can be observed for 1 — T and B^. 

7 Summary and conclusions 

In this paper we present determinations of the strong coupling as using event shape 
observable distributions at cm. energies between 91 and 209 GeV. Fits using NNLO and 
combined NNLO+NLLA predictions are used to extract as- Combining the results from 
the NNLO fits to the six event shape observables 1 — T, Mh, -£>w, -E>t, C and y^ 3 at the 

4 Thc terms inducing the renormalization scale dependence resemble terms of higher order in ag 
weighted with lna;^. 
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thirteen OPAL energy points yields 

a s (m z o) = 0.1201 ± 0.0008(stat.) ± 0.0013(exp.) ± 0.0010(had.) ± 0.0024(theo.) , 

with a total uncertainty on as(^z°) of 2.6%. Combining the results from the NNLO+NLLA 
fits yields 

a s (m z0 ) = 0.1189 ± 0.0008(stat.) ± 0.0016(exp.) ± 0.0010(had.) ± 0.0036(theo.) , 

with a total uncertainty of 3.4%. This result supersedes that presented in [TS] because 
it is based on more complete theoretical predictions while using the same experimental 
data and procedures. The variations between the combined results for as at the different 
energies are consistent with the running of as as predicted by QCD and exclude the 
absence of running. The as results we find for all 13 energy points are given in Tabs. [TDT - 
[13 for the NNLO, NNLO+NLLA, and NLO+NLLA calculations. 

The investigation of the renormalization scale dependence of as{m z o) shows a reduced 
dependence on when NNLO or NNLO+NLLA predictions are used, compared to anal- 
yses with NLO or NLO+NLLA predictions. The more complete NNLO or NNLO+NLLA 
QCD predictions thus lead to smaller theoretical uncertainties in our analysis. Adding 
the NLLA terms to the NNLO predictions significantly improves the description of the 
data. The standard procedure to quantify theoretical uncertainties leads, however, to 
somewhat larger uncertainties in most of the NNLO+NLLA fits. 
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Table 1: Year of data collection, energy range, mean cm. energy, integrated luminosity 
L and numbers of selected events for each OPAL data sample used in this analysis, see 
also [15] . The horizontal lines divide the data into four energy ranges used for presentation 
purposes. 



Year 




Range 


of y/s 


Mean A 




Selected 






[GeV] 




[GeV] 


[pb- 1 ] 


events 


1996, 


2000 


91.0- 


-91.5 


91.3 


14.7 


395695 


1995, 


1997 


129.9- 


-130.2 


130.1 


5.31 


318 


1995, 


1997 


136.0- 


-136.3 


136.1 


5.95 


312 


1996 




161.2- 


-161.6 


161.3 


10.06 


281 


1996 




170.2- 


-172.5 


172.1 


10.38 


218 


1997 




180.8- 


-184.2 


182.7 


57.72 


1077 


1998 




188.3- 


-189.1 


188.6 


185.2 


3086 


1999 




191.4- 


-192.1 


191.6 


29.53 


514 


1999 




195.4- 


-196.1 


195.5 


76.67 


1137 


1999, 2000 


199.1- 


-200.2 


199.5 


79.27 


1090 


1999, 


2000 


201.3- 


-202.1 


201.6 


37.75 


519 


2000 




202.5- 


-205.5 


204.9 


82.01 


1130 


2000 




205.5- 


-208.9 


206.6 


138.8 


1717 



Table 2: Fit ranges at all cm. energies. 



1 — T 


M H 




0.05-0.30 


0.17-0.45 


0.075-0.25 




C 




0.05-0.20 


0.18-0.60 


0.012-0.13 
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Table 3: Measurements of as using NNLO predictions and event shape distributions in 
four ranges of c.m. energy: at 91.3 GeV, 133.1 GeV, 161-183 GeV (177.4 GeV on average) 
and 189-209 GeV (197.0 GeV on average). 



i/i [GeV] Obs. «s(v^) ±stat. ±exp. ±had. ±theo. ±tot. 

91.3 1-T 0.1220 0.0002 0.0011 0.0030 0.0042 0.0053 

91.3 M H 0.1228 0.0002 0.0008 0.0026 0.0028 0.0039 

91.3 B T 0.1193 0.0002 0.0007 0.0033 0.0039 0.0052 

91.3 B w 0.1201 0.0001 0.0014 0.0010 0.0021 0.0027 

91.3 C 0.1188 0.0002 0.0009 0.0032 0.0035 0.0048 

91.3 «?, 0.1202 0.0002 0.0025 0.0005 0.0019 0.0032 



133.1 1-T 0.1126 0.0043 0.0038 0.0026 0.0032 0.0071 

133.1 M H 0.1110 0.0043 0.0033 0.0003 0.0020 0.0058 

133.1 B T 0.1065 0.0038 0.0048 0.0024 0.0027 0.0071 

133.1 B w 0.1123 0.0040 0.0026 0.0010 0.0017 0.0052 

133.1 C 0.1051 0.0046 0.0029 0.0033 0.0024 0.0068 

133.1 yg 0.1071 0.0056 0.0053 0.0014 0.0012 0.0079 

177.4 1-T 0.1088 0.0029 0.0027 0.0013 0.0028 0.0050 

177.4 M H 0.1081 0.0028 0.0042 0.0011 0.0018 0.0055 

177.4 B T 0.1051 0.0023 0.0031 0.0024 0.0026 0.0052 

177.4 B w 0.1047 0.0024 0.0029 0.0013 0.0012 0.0042 

177.4 C 0.1067 0.0028 0.0030 0.0015 0.0025 0.0050 

177.4 «?, 0.1084 0.0041 0.0031 0.0005 0.0013 0.0053 



197.0 1-T 0.1109 0.0012 0.0019 0.0011 0.0030 0.0039 

197.0 M H 0.1075 0.0012 0.0021 0.0016 0.0018 0.0034 

197.0 B T 0.1092 0.0011 0.0019 0.0013 0.0029 0.0039 

197.0 B w 0.1069 0.0011 0.0011 0.0004 0.0013 0.0021 

197.0 C 0.1086 0.0013 0.0017 0.0016 0.0026 0.0037 

197.0 t/P, 0.1073 0.0019 0.0023 0.0003 0.0013 0.0033 
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Table 4: Measurements of as using NNLO predictions and event shape distributions 
at 91.3 GeV. The employed covariance matrix includes an estimate of the experimental 
systematic covariance. The statistical error is taken from the fit with only the statisti- 
cal covariance. The experimental systematic uncertainty is given by the fit error with 
quadratically subtracted statistical uncertainty. 



Obs. 


a s (91.3 GeV) 


±stat. 


±exp. 


X 2 /d.o.f. 


l-T 


0.1168 


0.0002 


0.0007 


111/5 


M H 


0.1291 


0.0002 


0.0027 


18.3/4 


Bf 


0.1247 


0.0002 


0.0011 


311/5 


-Bw 


0.1107 


0.0001 


0.0009 


142/4 


C 


0.1177 


0.0002 


0.0012 


45.0/4 




0.1183 


0.0002 


0.0010 


1.6/3 
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Table 5: Measurements of as using NNLO+NLLA predictions and event shape distribu- 
tions in four ranges of cm. energy: at 91.3 GeV, 133.1 GeV, 161-183 GeV (177.4 GeV 
on average) and 189-209 GeV (197.0 GeV on average). 



i/i [GeV] Obs. «s(a/s) ±stat. ±exp. ±had. ±theo. ±tot. 

91.3 1-T 0.1219 0.0002 0.0012 0.0030 0.0041 0.0052 

91.3 M H 0.1207 0.0002 0.0008 0.0022 0.0033 0.0041 

91.3 B T 0.1213 0.0002 0.0010 0.0023 0.0048 0.0054 

91.3 B w 0.1164 0.0001 0.0013 0.0011 0.0041 0.0044 

91.3 C 0.1186 0.0002 0.0009 0.0030 0.0046 0.0056 

91.3 «?, 0.1195 0.0002 0.0025 0.0004 0.0023 0.0034 



133.1 1-T 0.1128 0.0044 0.0040 0.0025 0.0032 0.0072 

133.1 M H 0.1094 0.0042 0.0029 0.0003 0.0023 0.0056 

133.1 B T 0.1085 0.0041 0.0055 0.0016 0.0034 0.0078 

133.1 B w 0.1082 0.0035 0.0028 0.0013 0.0033 0.0057 

133.1 C 0.1049 0.0047 0.0034 0.0027 0.0030 0.0071 

133.1 yg 0.1068 0.0056 0.0051 0.0012 0.0017 0.0079 

177.4 1-T 0.1089 0.0029 0.0028 0.0012 0.0027 0.0050 

177.4 M H 0.1069 0.0027 0.0039 0.0013 0.0021 0.0053 

177.4 B T 0.1070 0.0026 0.0034 0.0019 0.0032 0.0057 

177.4 B w 0.1018 0.0023 0.0026 0.0012 0.0027 0.0046 

177.4 C 0.1063 0.0029 0.0031 0.0014 0.0030 0.0054 

177.4 «?, 0.1079 0.0040 0.0031 0.0005 0.0016 0.0053 



197.0 1-T 0.1110 0.0013 0.0020 0.0012 0.0029 0.0039 

197.0 M H 0.1064 0.0012 0.0021 0.0017 0.0021 0.0036 

197.0 B T 0.1112 0.0012 0.0021 0.0011 0.0035 0.0044 

197.0 B w 0.1045 0.0010 0.0013 0.0004 0.0029 0.0034 

197.0 C 0.1083 0.0013 0.0019 0.0017 0.0032 0.0043 

197.0 yPo 0.1068 0.0019 0.0022 0.0003 0.0016 0.0033 
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Table 6: Measurements of as using NNLO+NLLA predictions and event shape distribu- 
tions at 91.3 GeV. The employed covariance matrix includes an estimate of the exper- 
imental systematic covariance. The statistical error is taken from the fit with only the 
statistical covariance. The experimental systematic uncertainty is given by the fit error 
with quadratically subtracted statistical uncertainty. 



Obs. 


a s (91.3 GeV) 


±stat. 


±exp. 


X 2 /d.o.f. 


l-T 


0.1216 


0.0002 


0.0008 


17.0/5 


M H 


0.1215 


0.0002 


0.0021 


2.6/4 


Bf 


0.1224 


0.0002 


0.0009 


6.1/5 


-Bw 


0.1145 


0.0001 


0.0010 


3.4/4 


C 


0.1183 


0.0002 


0.0012 


4.7/4 




0.1191 


0.0002 


0.0011 


0.3/3 



Table 7: Combined values of ag(\/i) at the OPAL cm. energy ranges from NNLO (upper 
section) and NNLO+NLLA (lower section) analyses together with the statistical, experi- 
mental, hadronisation, theory and total uncertainties. 



[GeV] 


ats(y/s) 


±stat. 


±exp. 


±had. 


±theo. 


±tot. 


NNLO 


91.3 


0.1205 


0.0001 


0.0011 


0.0016 


0.0027 


0.0033 


133.1 


0.1109 


0.0036 


0.0027 


0.0011 


0.0022 


0.0051 


177.4 


0.1057 


0.0022 


0.0027 


0.0013 


0.0020 


0.0042 


197.0 


0.1077 


0.0010 


0.0012 


0.0006 


0.0018 


0.0025 


NNLO+NLLA 


91.3 


0.1196 


0.0002 


0.0012 


0.0013 


0.0036 


0.0040 


133.1 


0.1088 


0.0036 


0.0029 


0.0009 


0.0030 


0.0056 


177.4 


0.1042 


0.0024 


0.0027 


0.0009 


0.0027 


0.0046 


197.0 


0.1069 


0.0010 


0.0016 


0.0008 


0.0029 


0.0036 
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Table 8: Combined values of ag(mzo) at OPAL cm. energy ranges from NNLO (upper 
section) and NNLO+NLLA (lower section) analyses together with the statistical, experi- 
mental, hadronisation, theory and total uncertainties. 



v^[GeV] 


as(m z o) 


±stat. 


±exp. 


±had. 


±theo. 


±tot. 


NNLO 


130.1—206.6 
91.3—206.6 


0.1200 
0.1201 


0.0012 
0.0008 


0.0016 
0.0013 


0.0008 
0.0010 


0.0023 
0.0024 


0.0032 
0.0031 






NNLO- 


fNLLA 








130.1—206.6 
91.3—206.6 


0.1186 
0.1189 


0.0011 
0.0008 


0.0020 
0.0016 


0.0009 
0.0010 


0.0036 
0.0036 


0.0040 
0.0041 



Table 9: Combined values of a^mip) for each observable from NNLO (upper section) 
and NNLO+NLLA (lower section) analyses together with the statistical, experimental, 
hadronisation, theory and total uncertainties. 



Obs. 


a s (m z o) 


±stat. 


±exp. 


±had. 


±theo. 


±tot. 


NNLO 


1 — T 


0.1230 


0.0010 


0.0019 


0.0018 


0.0039 


0.0048 


M H 


0.1212 


0.0009 


0.0017 


0.0020 


0.0025 


0.0037 




0.1205 


0.0009 


0.0018 


0.0016 


0.0037 


0.0045 




0.1195 


0.0009 


0.0012 


0.0007 


0.0018 


0.0024 


C 


0.1199 


0.0011 


0.0016 


0.0023 


0.0033 


0.0045 




0.1199 


0.0008 


0.0026 


0.0004 


0.0018 


0.0033 


NNLO+NLLA 


1 — T 


0.1230 


0.0010 


0.0020 


0.0020 


0.0038 


0.0048 


M H 


0.1194 


0.0009 


0.0017 


0.0020 


0.0029 


0.0040 




0.1227 


0.0009 


0.0020 


0.0016 


0.0045 


0.0053 




0.1159 


0.0009 


0.0015 


0.0007 


0.0037 


0.0042 


C 


0.1196 


0.0011 


0.0019 


0.0023 


0.0042 


0.0053 




0.1192 


0.0010 


0.0026 


0.0004 


0.0021 


0.0035 
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Figure 1: Comparison of NNLO+NLLA calculations for the 1 — T, Mh, B^, B^, C and 
y® 3 variables with the parton level predictions of the Monte Carlo generators PYTHIA 6.1, 
HERWIG 6.2 and ARIADNE 4.11 at ^=91 GeV (see Sect. E3>- The vertically sym- 
metrical band between solid lines shows the maximum deviation of the ratio from one, 
between the three generators MQ in the positive and the negative direction. The arrows 
indicate the respective fit ranges. 
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Figure 2: The points show the 1 — T and Mh distributions at the hadron level for y/s = 
91.3 and 206.6 GeV with statistical uncertainty bars. Some uncertainty bars are smaller 
than the data points. Superimposed as histograms are the NNLO and NNLO+NLLA 
predictions combined with hadronisation effects using the corresponding fit results for 
&s(\/s) shown in Tabs. El |5j The arrows indicate the fit ranges. 




Figure 3: The points show the and B-w distributions at the hadron level for y/s = 
91.3 and 206.6 GeV with statistical uncertainty bars. Some uncertainty bars are smaller 
than the data points. Superimposed as histograms are the NNLO and NNLO+NLLA 
predictions combined with hadronisation effects using the corresponding fit results for 
ocs(y/s) shown in Tabs. El |5j The arrows indicate the fit ranges. 
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Figure 4: The points show the C and distributions at the hadron level for yfs = 
91.3 and 206.6 GeV with statistical uncertainty bars. Some uncertainty bars are smaller 
than the data points. Superimposed as histograms are the NNLO and NNLO+NLLA 
predictions combined with hadronisation effects using the corresponding fit results for 
a s(Vs) shown in Tabs. El El The arrows indicate the fit ranges. 
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25 50 75 100 125 150 175 200 

Vs[GeV] 

Figure 5: The points show the values of as for the OPAL energy ranges. The inner 
uncertainty bars show the combined statistical and experimental uncertainties and the 
outer the total uncertainties. The full and dashed lines indicate the as result from the 
NNLO+NLLA analysis that combines all variables and OPAL energy points. The results 
from the NNLO+NLLA analysis of JADE data [12] are shown as well. 
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Figure 6: as results combined over all OPAL cm. energies for different event shape vari- 
ables and different QCD calculations as indicated on the figure. The shaded bands and 
dashed lines show the values of as(mzo) combined from these values with total uncertain- 
ties. The inner and outer uncertainty bars show the combined statistical and experimental 
and total uncertainties, respectively. 
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Figure 7: Results for as and x 2 /d.o.f. for fits of QCD predictions to the C distribution, 
as a function of x„. 
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A Appendix: Detailed results 



Table 10: Results of NNLO fits to event shape observable distributions at the OPAL cm. 
energies. The x 2 /d.o.f. values are based on the statistical errors only. 



a/s [GeV] Obs. as(Vs) istat. ±exp. ±had. ±theo. x 2 /d-o.f. 

91.3 1-T 0.1220 0.0002 0.0011 0.0030 0.0042 440.1/5 

91.3 M H 0.1228 0.0002 0.0008 0.0027 0.0028 393.8/4 

91.3 B T 0.1193 0.0001 0.0007 0.0033 0.0039 1142.4/5 

91.3 B w 0.1201 0.0001 0.0013 0.0010 0.0021 446.6/4 

91.3 C 0.1188 0.0002 0.0009 0.0032 0.0035 531.8/4 

91.3 yg 0.1202 0.0002 0.0025 0.0005 0.0019 31.5/3 

130.1 1-T 0.1179 0.0057 0.0030 0.0025 0.0037 8.3/5 

130.1 M H 0.1158 0.0057 0.0027 0.0004 0.0023 6.5/4 

130.1 B T 0.1098 0.0051 0.0045 0.0019 0.0029 8.2/5 

130.1 B w 0.1155 0.0049 0.0020 0.0006 0.0018 0.8/4 

130.1 C 0.1106 0.0061 0.0019 0.0028 0.0027 7.8/4 

130.1 yg 0.1106 0.0075 0.0046 0.0016 0.0014 4.7/3 

136.1 1-T 0.1019 0.0062 0.0080 0.0029 0.0022 12.0/5 

136.1 M H 0.1012 0.0061 0.0067 0.0016 0.0014 5.1/4 

136.1 B T 0.1002 0.0053 0.0072 0.0034 0.0023 8.7/5 

136.1 B w 0.1021 0.0051 0.0082 0.0024 0.0013 7.0/4 

136.1 C 0.0932 0.0065 0.0073 0.0045 0.0016 10.6/4 

136.1 yg 0.0999 0.0076 0.0095 0.0009 0.0010 3.0/3 

161.3 1-T 0.1068 0.0065 0.0031 0.0025 0.0026 6.4/5 

161.3 M H 0.1073 0.0064 0.0040 0.0009 0.0017 3.9/4 

161.3 B T 0.1001 0.0055 0.0040 0.0034 0.0022 8.7/5 

161.3 B w 0.1022 0.0054 0.0023 0.0034 0.0012 9.4/4 

161.3 C 0.1019 0.0066 0.0044 0.0034 0.0021 8.2/4 

161.3 yg 0.1060 0.0086 0.0020 0.0008 0.0013 2.0/3 

172.1 1-T 0.1087 0.0076 0.0072 0.0030 0.0029 6.3/5 

172.1 M H 0.1068 0.0073 0.0062 0.0009 0.0018 5.0/4 

172.1 B T 0.1018 0.0064 0.0047 0.0021 0.0023 4.5/5 

172.1 B w 0.1001 0.0064 0.0039 0.0007 0.0011 1.9/4 

172.1 C 0.1046 0.0075 0.0071 0.0033 0.0024 7.0/4 

172.1 yg 0.1044 0.0104 0.0132 0.0013 0.0012 2.6/3 

182.7 T^T 0.1098 0.0035 0.0025 0.0018 O0029 0.6/5 

182.7 M H 0.1091 0.0033 0.0044 0.0015 0.0018 4.4/4 

182.7 B T 0.1078 0.0028 0.0029 0.0021 0.0028 9.0/5 

182.7 B w 0.1071 0.0029 0.0032 0.0005 0.0013 4.3/4 

182.7 C 0.1086 0.0034 0.0026 0.0018 0.0026 1.1/4 

182.7 yg 0.1103 0.0048 0.0040 0.0003 0.0014 1.4/3 
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Table 11: Results of NNLO fits to event shape observable distributions at the OPAL cm. 
energies. The x 2 /d.o.f. values are based on the statistical errors only. 



t/s [GeV] Obs. tts{y/s) dbstat. zbexp. zbhad. zbtheo. x 2 /d.o.f. 

188.6 1-T 0.1133 0.0020 0.0015 0.0014 0.0032 14.3/5 

188.6 M H 0.1090 0.0020 0.0021 0.0018 0.0018 4.9/4 

188.6 B T 0.1100 0.0017 0.0013 0.0011 0.0030 5.0/5 

188.6 B w 0.1074 0.0017 0.0006 0.0006 0.0014 2.1/4 

188.6 C 0.1094 0.0021 0.0012 0.0016 0.0027 6.4/4 

188.6 yg 0.1083 0.0029 0.0020 0.0004 0.0013 3.5/3 

191.6 1-T 0.1112 0.0050 0.0085 0.0017 0.0031 1.9/5 

191.6 M H 0.1094 0.0048 0.0055 0.0014 0.0019 4.1/4 

191.6 B T 0.1015 0.0046 0.0089 0.0025 0.0022 10.2/5 

191.6 B w 0.1044 0.0042 0.0050 0.0005 0.0013 0.8/4 

191.6 C 0.1053 0.0053 0.0060 0.0020 0.0024 0.3/4 

191.6 yg 0.1051 0.0071 0.0066 0.0005 0.0012 0.0/3 

195.5 1-T 0.1119 0.0033 0.0040 0.0021 0.0031 9.6/5 

195.5 A/jj 0.1057 0.0033 0.0032 0.0018 0.0016 2.0/4 

195.5 B T 0.1081 0.0029 0.0026 0.0016 0.0028 6.7/5 

195.5 B w 0.1062 0.0029 0.0018 0.0007 0.0013 4.5/4 

195.5 C 0.1087 0.0035 0.0035 0.0017 0.0026 1.9/4 

195.5 yg 0.1008 0.0048 0.0046 0.0004 0.0010 1.4/3 

199.5 1-T 0.1105 0.0034 0.0057 0.0021 0.0029 8.5/5 

199.5 M H 0.1061 0.0034 0.0028 0.0007 0.0017 7.1/4 

199.5 B T 0.1102 0.0030 0.0039 0.0023 0.0030 8.5/5 

199.5 B w 0.1076 0.0029 0.0020 0.0010 0.0013 3.9/4 

199.5 C 0.1077 0.0036 0.0039 0.0012 0.0026 5.5/4 

199.5 yg 0.1130 0.0049 0.0045 0.0009 0.0015 4.8/3 

201.6 T^T 0.1168 0.0048 0.0051 0.0013 01)035 1.5/5 
201.6 M H 0.1045 0.0049 0.0047 0.0012 0.0016 12.2/4 
201.6 B T 0.1072 0.0044 0.0027 0.0015 0.0027 8.4/5 
201.6 B w 0.1057 0.0043 0.0026 0.0017 0.0012 7.7/4 
201.6 C 0.1092 0.0052 0.0061 0.0024 0.0026 11.3/4 
201.6 yg 0.1107 0.0071 0.0076 0.0008 0.0014 1.4/3 
204.9 1 -T 0.1111 0.0032 0.0032 0.0015 0.0031 8.7/5 
204.9 M H 0.1134 0.0029 0.0036 0.0006 0.0021 7.0/4 
204.9 B T 0.1091 0.0028 0.0057 0.0011 0.0029 4.7/5 
204.9 B w 0.1108 0.0027 0.0034 0.0005 0.0015 3.1/4 
204.9 C 0.1085 0.0033 0.0068 0.0017 0.0027 3.3/4 
204.9 yg 0.1090 0.0045 0.0060 0.0007 0.0013 3.7/3 
206.6 1 -T 0.1066 0.0026 0.0019 0.0018 0.0026 2.9/5 
206.6 M H 0.1055 0.0025 0.0017 0.0021 0.0016 0.8/4 
206.6 B T 0.1100 0.0021 0.0027 0.0010 0.0030 0.8/5 
206.6 B w 0.1061 0.0022 0.0024 0.0006 0.0013 0.9/4 
206.6 C 0.1079 0.0026 0.0020 0.0017 0.0026 1.6/4 
206.6 yg 0.1058 0.0037 0.0026 0.0003 0.0012 0.6/3 
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Table 12: Results of NNLO+NLLA fits to event shape observable distributions at the 
OPAL cm. energies. The x 2 /d.o.i. values are based on the statistical errors only. 



y/s [GeV] Obs. as(v^) ±stat. ±exp. zbhad. zbtheo. x 2 /d.o.f. 

91.3 1-T 0.1219 0.0002 0.0012 0.0030 0.0041 120.5/5 

91.3 M B 0.1206 0.0002 0.0008 0.0022 0.0033 66.8/4 

91.3 B T 0.1213 0.0002 0.0009 0.0023 0.0048 40.0/5 

91.3 B w 0.1164 0.0001 0.0013 0.0011 0.0041 108.9/4 

91.3 C 0.1187 0.0002 0.0009 0.0030 0.0046 51.2/4 

91.3 yg 0.1195 0.0002 0.0025 0.0005 0.0023 8.0/3 

130.1 1-T 0.1178 0.0058 0.0031 0.0025 0.0037 8.1/5 

130.1 M H 0.1138 0.0054 0.0022 0.0008 0.0026 6.7/4 

130.1 B T 0.1130 0.0056 0.0050 0.0024 0.0039 5.5/5 

130.1 B w 0.1119 0.0046 0.0019 0.0009 0.0036 1.0/4 

130.1 C 0.1113 0.0063 0.0024 0.0019 0.0038 6.4/4 

130.1 yg 0.1105 0.0075 0.0044 0.0015 0.0019 4.4/3 

136.1 1-T 0.1026 0.0063 0.0084 0.0025 0.0023 11.0/5 

136.1 M H 0.1003 0.0059 0.0066 0.0018 0.0017 4.7/4 

136.1 B T 0.1012 0.0056 0.0079 0.0034 0.0026 8.2/5 

136.1 B w 0.0993 0.0047 0.0076 0.0024 0.0024 7.3/4 

136.1 C 0.0927 0.0064 0.0073 0.0044 0.0016 10.5/4 

136.1 yg 0.0996 0.0075 0.0093 0.0008 0.0012 2.9/3 

161.3 1-T 0.1076 0.0066 0.0032 0.0020 0.0028 5.6/5 

161.3 M H 0.1068 0.0062 0.0037 0.0017 0.0022 3.1/4 

161.3 B T 0.1021 0.0059 0.0041 0.0029 0.0027 7.5/5 

161.3 B w 0.1001 0.0051 0.0022 0.0030 0.0025 9.1/4 

161.3 C 0.1024 0.0067 0.0046 0.0029 0.0028 7.4/4 

161.3 yg 0.1057 0.0085 0.0020 0.0007 0.0015 1.9/3 

172.1 1-T 0.1080 0.0076 0.0071 0.0032 0.0024 6.7/5 

172.1 M H 0.1047 0.0069 0.0057 0.0015 0.0019 5.6/4 

172.1 B T 0.1038 0.0069 0.0054 0.0015 0.0029 3.7/5 

172.1 B w 0.0978 0.0060 0.0038 0.0007 0.0023 2.0/4 

172.1 C 0.1038 0.0075 0.0072 0.0034 0.0025 7.2/4 

172.1 yg 0.1038 0.0102 0.0129 0.0013 0.0015 2.6/3 

182.7 1-T 0.1098 0.0035 0.0026 0.0018 0.0028 0.4/5 

182.7 M H 0.1080 0.0032 0.0040 0.0017 0.0022 3.7/4 

182.7 B T 0.1102 0.0031 0.0032 0.0015 0.0035 5.5/5 

182.7 B w 0.1047 0.0027 0.0030 0.0004 0.0029 3.5/4 

182.7 C 0.1083 0.0035 0.0028 0.0019 0.0032 0.6/4 

182.7 yg 0.1098 0.0047 0.0039 0.0003 0.0017 1.3/3 
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Table 13: Results of NNLO+NLLA fits to event shape observable distributions at the 
OPAL cm. energies. The x 2 /d.o.f. values are based on the statistical errors only. 



t/s [GeV] Obs. cts{y/s) ±stat. zbexp. zbhad. zbtheo. x 2 /d.o.f. 

188.6 1-T 0.1133 0.0020 0.0016 0.0012 0.0032 13.6/5 

188.6 M H 0.1077 0.0019 0.0020 0.0019 0.0021 4.3/4 

188.6 B T 0.1119 0.0019 0.0014 0.0012 0.0036 2.1/5 

188.6 B w 0.1048 0.0016 0.0006 0.0006 0.0029 1.1/4 

188.6 C 0.1092 0.0021 0.0014 0.0018 0.0033 4.0/4 

188.6 yg 0.1077 0.0028 0.0020 0.0004 0.0016 3.6/3 

191.6 1-T 0.1110 0.0051 0.0085 0.0016 0.0028 2.1/5 

191.6 M H 0.1076 0.0046 0.0055 0.0011 0.0021 4.8/4 

191.6 B T 0.1040 0.0050 0.0099 0.0021 0.0028 7.9/5 

191.6 B w 0.1019 0.0040 0.0049 0.0005 0.0026 0.9/4 

191.6 C 0.1047 0.0053 0.0060 0.0020 0.0027 0.3/4 

191.6 yg 0.1046 0.0070 0.0066 0.0005 0.0014 0.0/3 

195.5 1-T 0.1121 0.0034 0.0041 0.0020 0.0031 9.0/5 

195.5 M H 0.1047 0.0032 0.0032 0.0020 0.0020 1.2/4 

195.5 B T 0.1104 0.0032 0.0028 0.0012 0.0035 3.9/5 

195.5 B w 0.1038 0.0027 0.0017 0.0005 0.0028 3.7/4 

195.5 C 0.1083 0.0036 0.0034 0.0018 0.0031 1.6/4 

195.5 yg 0.1005 0.0047 0.0046 0.0003 0.0013 1.2/3 

199.5 1-T 0.1109 0.0035 0.0059 0.0017 0.0031 7.2/5 

199.5 M H 0.1052 0.0033 0.0029 0.0010 0.0020 6.0/4 

199.5 B T 0.1129 0.0034 0.0045 0.0019 0.0039 4.9/5 

199.5 B w 0.1053 0.0028 0.0020 0.0006 0.0029 2.5/4 

199.5 C 0.1077 0.0037 0.0040 0.0013 0.0032 4.2/4 

199.5 yg 0.1124 0.0048 0.0044 0.0009 0.0019 4.9/3 

201.6 T^T 0.1172 0.0049 0.0053 0.0010 (X0037 1.0/5 
201.6 M H 0.1036 0.0048 0.0047 0.0010 0.0019 11.7/4 
201.6 B T 0.1097 0.0048 0.0036 0.0011 0.0034 6.5/5 
201.6 B w 0.1035 0.0041 0.0025 0.0013 0.0028 7.1/4 
201.6 C 0.1104 0.0054 0.0068 0.0020 0.0038 9.3/4 
201.6 yg 0.1104 0.0070 0.0074 0.0007 0.0018 1.2/3 
204.9 1 -T 0.1108 0.0032 0.0031 0.0014 0.0028 9.9/5 
204.9 M H 0.1117 0.0029 0.0034 0.0005 0.0024 7.9/4 
204.9 B T 0.1108 0.0030 0.0060 0.0011 0.0034 4.5/5 
204.9 B w 0.1077 0.0025 0.0032 0.0004 0.0032 4.7/4 
204.9 C 0.1081 0.0033 0.0070 0.0017 0.0031 3.3/4 
204.9 yg 0.1085 0.0044 0.0060 0.0007 0.0017 3.5/3 
206.6 1 -T 0.1068 0.0027 0.0019 0.0021 0.0025 2.1/5 
206.6 M H 0.1042 0.0025 0.0017 0.0022 0.0019 0.7/4 
206.6 B T 0.1116 0.0023 0.0030 0.0012 0.0035 1.5/5 
206.6 B w 0.1035 0.0021 0.0022 0.0005 0.0028 1.4/4 
206.6 C 0.1076 0.0027 0.0021 0.0018 0.0031 0.9/4 
206.6 yg 0.1054 0.0036 0.0025 0.0004 0.0015 0.5/3 
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Table 14: Results of NLO+NLLA fits to event shape observable distributions at the OPAL 
cm. energies. The x 2 /d.o.i. values are based on the statistical errors only. 



y/s [GeV] Obs. asly/s) ±stat. zbexp. zbhad. zbtheo. x 2 / d -o.f. 

91.3 1-T 0.1233 0.0002 0.0012 0.0031 0.0073 247.1/5 

91.3 M B 0.1196 0.0002 0.0008 0.0020 0.0057 36.9/4 

91.3 B T 0.1226 0.0002 0.0009 0.0029 0.0076 190.5/5 

91.3 B w 0.1145 0.0001 0.0013 0.0012 0.0055 140.2/4 

91.3 C 0.1173 0.0002 0.0009 0.0031 0.0064 144.3/4 

91.3 yg 0.1192 0.0002 0.0025 0.0005 0.0044 4.5/3 

130.1 1-T 0.1187 0.0059 0.0032 0.0027 0.0065 8.6/5 

130.1 M H 0.1129 0.0054 0.0021 0.0008 0.0048 6.4/4 

130.1 B T 0.1138 0.0057 0.0050 0.0024 0.0063 5.8/5 

130.1 B w 0.1100 0.0046 0.0019 0.0010 0.0050 1.1/4 

130.1 C 0.1100 0.0061 0.0023 0.0022 0.0054 6.8/4 

130.1 yg 0.1104 0.0074 0.0043 0.0014 0.0035 4.2/3 

136.1 1-T 0.1031 0.0064 0.0085 0.0027 0.0044 11.6/5 

136.1 M H 0.0995 0.0058 0.0066 0.0018 0.0032 4.4/4 

136.1 B T 0.1017 0.0057 0.0079 0.0036 0.0048 8.7/5 

136.1 B w 0.0976 0.0047 0.0076 0.0024 0.0036 7.3/4 

136.1 C 0.0919 0.0063 0.0071 0.0044 0.0030 10.6/4 

136.1 yg 0.0995 0.0075 0.0094 0.0008 0.0026 2.9/3 

161.3 1-T 0.1083 0.0067 0.0032 0.0023 0.0051 5.9/5 

161.3 M H 0.1061 0.0062 0.0035 0.0019 0.0041 2.7/4 

161.3 B T 0.1026 0.0060 0.0041 0.0032 0.0049 7.8/5 

161.3 B w 0.0984 0.0051 0.0022 0.0029 0.0037 8.9/4 

161.3 C 0.1012 0.0065 0.0045 0.0030 0.0042 7.7/4 

161.3 yg 0.1056 0.0085 0.0019 0.0006 0.0030 1.8/3 

172.1 1-T 0.1090 0.0077 0.0073 0.0033 0.0050 6.6/5 

172.1 M H 0.1036 0.0069 0.0056 0.0016 0.0035 5.8/4 

172.1 B T 0.1045 0.0070 0.0054 0.0016 0.0051 3.7/5 

172.1 B w 0.0961 0.0060 0.0037 0.0007 0.0035 2.0/4 

172.1 C 0.1029 0.0074 0.0070 0.0034 0.0042 7.2/4 

172.1 yg 0.1038 0.0102 0.0130 0.0012 0.0029 2.6/3 

182.7 1-T 0.1108 0.0036 0.0026 0.0018 0.0053 0.6/5 

182.7 M H 0.1071 0.0032 0.0039 0.0018 0.0041 3.6/4 

182.7 B T 0.1110 0.0031 0.0033 0.0017 0.0059 6.0/5 

182.7 B w 0.1029 0.0027 0.0029 0.0004 0.0042 3.6/4 

182.7 C 0.1073 0.0034 0.0027 0.0019 0.0049 0.6/4 

182.7 yg 0.1096 0.0047 0.0039 0.0003 0.0034 1.3/3 
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Table 15: Results of NLO+NLLA fits to event shape observable distributions at the OPAL 
cm. energies. The x 2 /d.o.i. values are based on the statistical errors only. 



t/s [GeV] Obs. tts{y/s) dbstat. zbexp. zbhad. zbtheo. x 2 /d.o.f. 

188.6 l-T 0.1143 0.0021 0.0016 0.0014 0.0058 14.2/5 

188.6 M H 0.1068 0.0019 0.0020 0.0019 0.0040 4.1/4 

188.6 B T 0.1128 0.0019 0.0014 0.0012 0.0061 2.7/5 

188.6 B w 0.1031 0.0016 0.0006 0.0006 0.0042 1.2/4 

188.6 C 0.1083 0.0021 0.0013 0.0018 0.0050 4.1/4 

188.6 yg 0.1075 0.0028 0.0020 0.0004 0.0032 3.4/3 

191.6 l-T 0.1120 0.0052 0.0087 0.0016 0.0055 2.0/5 

191.6 M H 0.1067 0.0046 0.0055 0.0010 0.0039 5.0/4 

191.6 B T 0.1045 0.0051 0.0100 0.0024 0.0051 8.4/5 

191.6 B w 0.1002 0.0040 0.0048 0.0004 0.0039 0.9/4 

191.6 C 0.1039 0.0053 0.0059 0.0020 0.0044 0.3/4 

191.6 yg 0.1044 0.0070 0.0066 0.0005 0.0030 0.0/3 

195.5 l-T 0.1130 0.0034 0.0041 0.0021 0.0057 9.3/5 

195.5 A/jj 0.1040 0.0032 0.0032 0.0021 0.0038 0.8/4 

195.5 B T 0.1112 0.0032 0.0028 0.0013 0.0059 4.5/5 

195.5 B w 0.1021 0.0027 0.0017 0.0005 0.0041 3.7/4 

195.5 C 0.1073 0.0035 0.0033 0.0018 0.0049 1.7/4 

195.5 yg 0.1003 0.0047 0.0046 0.0003 0.0026 1.2/3 

199.5 l-T 0.1117 0.0036 0.0059 0.0020 0.0055 7.8/5 

199.5 M H 0.1045 0.0033 0.0029 0.0012 0.0039 5.2/4 

199.5 B T 0.1136 0.0034 0.0046 0.0022 0.0062 6.0/5 

199.5 B w 0.1036 0.0028 0.0019 0.0005 0.0042 2.2/4 

199.5 C 0.1066 0.0036 0.0039 0.0013 0.0049 4.6/4 

199.5 yg 0.1121 0.0048 0.0043 0.0009 0.0037 4.8/3 

201.6 T^T 0.1182 0.0050 0.0054 0.0011 0l)065 1.1/5 
201.6 M H 0.1029 0.0048 0.0048 0.0009 0.0036 11.4/4 
201.6 B T 0.1105 0.0049 0.0034 0.0011 0.0058 6.7/5 
201.6 B w 0.1018 0.0040 0.0025 0.0013 0.0040 7.1/4 
201.6 C 0.1091 0.0053 0.0065 0.0020 0.0053 9.7/4 
201.6 yg 0.1102 0.0070 0.0073 0.0007 0.0035 1.2/3 
204.9 l-T 0.1118 0.0032 0.0032 0.0014 0.0054 9.6/5 
204.9 M H 0.1109 0.0029 0.0034 0.0006 0.0045 7.9/4 
204.9 B T 0.1117 0.0031 0.0061 0.0011 0.0060 4.2/5 
204.9 B w 0.1059 0.0025 0.0031 0.0005 0.0045 5.2/4 
204.9 C 0.1072 0.0032 0.0069 0.0017 0.0049 3.0/4 
204.9 yg 0.1084 0.0044 0.0060 0.0006 0.0033 3.4/3 
206.6 l-T 0.1076 0.0027 0.0019 0.0020 0.0049 2.2/5 
206.6 M H 0.1034 0.0025 0.0017 0.0022 0.0036 0.8/4 
206.6 B T 0.1126 0.0024 0.0030 0.0012 0.0061 1.2/5 
206.6 B w 0.1018 0.0021 0.0022 0.0005 0.0040 1.5/4 
206.6 C 0.1066 0.0026 0.0020 0.0018 0.0048 1.1/4 
206.6 yg 0.1052 0.0036 0.0025 0.0004 0.0030 0.5/3 
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